Method for monopole creation in spinor Bose-Einstein condensates 
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We present a method for creating a monopole in an antiferromagnetic spin-1 Bose-Einstein con- 
densate. The required phase engineering of the multicomponent condensate is achieved using light 
shifts, which depend on both the magnetic substate mp and polarization of the incident laser beam. 
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I. INTRODUCTION 

Since the experimental realization of Bose-Einstein 
condensates (BEC's) in 1995 |,|, |, a wide variety of 
properties and effects in BEC's has been studied. Among 
the major topics of interest are topological defects such as 
vortices [| |, 0, 1 1, |l3 and solitons 0, |l|, |l||l| . More 
recently, the realization of BEC's in optical traps |l^ 
has opened up interesting new avenues of research in this 
area. This development makes it possible to create mul- 
ticomponent condensates where all magnetic substates 
mp are trapped. The resulting spin degree of freedom 
permits the existence of topological defects that cannot 
appear in scalar condensates, such as skyrmions and 
monopoles IT^ . Important properties of monopoles, such 
as stability and dynamics, have already been examined 
p?! , 18 . Here we propose an explicit method for creat- 
ing a monopole in BEC's, focusing on the experimentally 
relevant case of ^'^Na. Our proposal extends the method 
of "phase imprinting" used to generate topological de- 
fects in scalar condensates [|[ ^ |l^. In particular, we 
demonstrate that the energy shift experienced by each 
component of the condensate depends on both mp and 
the polarization of the laser beam, and that by apply- 
ing a sequence of different polarizations the phases of the 
components can be engineered in a broad manner. 

A spin-1 BEC can be described by the condensate 
wavefunction 'Fi — y^n{r)Q{r), where n{r) is the total 
density and C(r) is a normalized spinor describing the 
local spin of the condensate. For an anti-ferromagnetic 
BEC such as ^■^Na, the ground state spinor is given by flfl] 



C.(r) 



(1) 



The spinor corresponding to a monopole has the 
form (l| 



-r7ij;(r) + iniyir) 
C™(r) = ^ I V2m.(r) 



V2 



(2) 



m^{r) + imy{v) 



m(r) is a unit vector field that describes the local ro- 
tation of the ground state spinor. By considering the 
energy and winding number of such a condensate, it can 
be shown |jl^ that m(r) = v/r uniquely characterizes 
the monopole structure. We see that m(r) assumes the 
familiar hedgehog structure associated with a monopole. 

For the purpose of phase engineering, it is convenient 
to also write this spinor in spherical coordinates as 



1 



sm We 



C™(r) = -= \/2cos0 
\ sin6'e*'^ 



(3) 



where (p is the azimuthal angle and 9 the polar an- 
gle. In addition, we define the phase of each component 
Lpi = arg(Ci). For example, the phases of the components 
in the monopole are given by 

/ -0 + 7r\ 
^™ = I ^0{-^z) I , (4) 

where 0{z) is the step function. 

II. CONDENSATE DYNAMICS 

The Hamiltonian of the condensate in second quan- 
tized form is [Rl, Eo| 



dr f ^Vt^J • \/4>, + 14rap(r)V'l^» + 



(5) 
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where is the field annihilation operator for spin 

state z, Fi {i — x,y,z} are the spin matrices, Vtrap(r) is 
the trapping potential, and Cq = Airh'^^aQ + 2a2)/3m and 
C2 = 4:Trh^{a2 — ao)/3rn are constants that depend on the 
scattering lengths oq and 02 . Ti.B and TYl describe respec- 
tively the applied magnetic field and the laser beams that 
will be used to create the monopole. During the times 
that the applied magnetic field and lasers are turned on, 
we assume that their respective Hamiltonians dominate 
Ti and are the only terms that need to be considered in 
the evolution of the system. 
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The Hamiltonian for the apphed magnetic field is given 
by Ti-B = — J dr^l{'j'B ■ hFij)ipj, where 7 is the gyro- 
magnetic ratio. When we ignore all other terms in Ti., the 
magnetic field leads to the following equations of motion: 



a -polarized beams are 



(6) 



When a laser beam is applied, the ground state energy 
of magnetic substate i is shifted by an amount pl| 



-I X 



E 



f<2 



(7) 



where the sum is carried out over all excited states 
Here, (5^ is the detuning of the laser from the \i) \j) 
transition, / the laser intensity, T the rate of sponta- 
neous decay, loq the atomic resonance frequency, and Cy 
is a transition strength characterizing the coupling be- 
tween atom and laser. The coefficients Cij depend on 
the ground substate mp and on the polarization of the 
laser beam and are given in Figure |l|. 



III. PHASE ENGINEERING USING PHASE 
IMPRINTING 



The concept of phase engineering via imprinting was 
suggested |9[ ^ and experimentally demonstrated 

as a means of creating topological defects in scalar 
condensates. We extend the method to spinor conden- 
sates, focusing on the case of ^^Na. We show that the 
energy shifts induced by different beam polarizations al- 
low for wide capabilities in phase engineering. 

Suppose that a laser is detuned by an amount 6 from 
the '^Si/2iF = l)-^^Pi/2(i^ = 2) transition, and that the 
states '^Pi/2{F — 1) and '^Pi/2{F = 2) are separated in 
frequency by an amount rjT, as shown in Figure |^. From 
Eq. §, the energy of the ground state is shifted by an 
amount 




^ . (8) 



AE^ = 5-J X 



We assume that the laser pulse has a square envelope 
with a duration t that is short compared to the conden- 
sate's correlation time = h/ fi, where fj, is the chemical 
potential. In this case the primary effect of the light shift 



is to imprint a phase onto each component |12 



(9) 



Using the coefficients given in Figure |i| and the ex- 
pansion = 3- ^ ^ + • • 'j one can obtain explicit 
results for the phase shifts generated in the condensate. 
For S ^ rjT, the shifts corresponding to tt-, ct"''-, and 



A^+{tf 
A^-{tf 



rjT 



-at[{2 2 2)-^(l 1)], (10) 



1), 
3). 



(11) 
(12) 



Here a = nc^TI /4hu)§d. Equations (|l0[)-(gj) suggest 
a basic technique for phase engineering of the spinor 
condensate. To shift the phase of m^? = 1 relative to 
mp = —1, one can use beams with polarizations a'^ 
and cr^. Once the phase between these two states is 
fixed as desired, one can use 7r-polarized light to adjust 
the phase of mi? = relative to the other states. Al- 
though it is not used in our proposal for monopole cre- 
ation, we note that this phase shift is due to a second 
order term in Eq. (p^), and we thus expect that it will 
be more difficult to implement in an experiment. It may 
be useful, in addition, to alter the global phase of the 
condensate, which can be realized approximately by us- 
ing the first order term in Eq. (p^. Finally, one can 
incorporate spatial dependence into the phases by using 
the above pulses in conjunction with absorption plates 
to create non- uniform intensity profiles /(r) Writing 
/(r) = /op(r), where Iq is the maximum intensity in the 
profile and < p(r) < 1, Eqs. (pO|)-(p^) remain valid with 
the replacements a ~^ aQp{r) and aQ = Trc^TlQ/AhujQS. 



IV. MODEL FOR MONOPOLE CREATION 

Martikainen et al. point out that the monopole 
spinor given in Eq. (|^) resembles a combination of other 
topological defects. Indeed, the form of Cm(r) presented 
in Eq. makes it clear that the mj? = +1 state con- 
tains a vortex, the mp = — 1 state an anti- vortex, and 
the mp — state a soliton. They thus suggest that 
one can experimentally create monopoles by creating the 
appropriate defects in the individual components. Here 
we provide a method to accomplish this using the phase 
imprinting concepts discussed above. 

Starting from the ground state of Eq. (^ , we first ap- 
ply a constant magnetic field of magnitude Bq along the 
negative y-axis for a time ts- Defining ujb = jBq, Eq. (|^) 
readily yields 



-Sm{uJBtB) 

Cb = C(*s) = ^ ( V2cos{u}BtB) 
sm{ujBtB) 



1 



V2 



(13) 



The phases of (b (for < i^BtB < 7i'/2) are 

Vl^in 0). (14) 

The effect of the magnetic field is to populate the mp — 
±1 states prior to phase engineering. Experimentally, 
tB (0 < (jJBtB < 7r/2) should be chosen to best reahze 
the monopole structure. 
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To create the monopole, we assume that we have bhie- 
detuned lasers, all of which have equal intensities /q and 
values of S, and which can be applied in square-shaped 
pulses. The results can be accordingly modified to cover 
more general situations. In addition, we assume that 
we have absorption plates that can create the following 
intensity profiles: 



P2(r) 
P4(r) 



27r' 
(27r - cj)) 

e{-z). 



(15) 

(16) 

(17) 
(18) 



The pulse sequence to create the monopole is given in 
Table |. Pulses A and B use intensity profiles with az- 
imuthal dependence to create a vortex in the mp = 1 
component and an anti- vortex in the mj? = — 1 compo- 
nent. Pulses C and D together imprint the necessary 
phases to generate the soliton in the mp = component. 
It can be readily verified that the resulting phases of the 
components are exactly those of the monopole given in 
Eq. (^). The total time of the laser pulses, assuming that 
pulses with profiles P3(r) and pi(j) can be executed si- 
multaneously, is T = Sn/aQ. Strictly speaking, pulse C 
is not physically possible — the profile (r) requires an 
axis of propagation k in the x^y plane, which is inconsis- 
tent with a polarization vector of e = — . However, it 
is possible to create an effective interaction term resem- 
bling e • r using a beam that does propagate in the x—y 
plane. Indeed, consider a pulse with polarization vec- 
tor e' = — sandwiched between magnetic field pulses 

that generate the rotations exp {—i^Fy) and exp f -Fj,) ■ 
One can readily check that 



exp 



r exp I I 



(19) 



Physically, the m p eigenstates along the z-axis are trans- 
formed by rotation to become mp eigenstates along the 
X-axis. The phase shifts are subsequently performed on 
these states, where P3,{y) and e' are consistent, before the 
system is rotated back. 

Experimentally, one should seek to make 5 as large as 
possible while increasing /q to maintain a given potential 
depth. This is because the photon scattering rate Psc of 
an individual atom scales as /o(r/(5)^ compared to I^^T / 5 
for the dipole potential pi| , and also because the approx- 
imation errors in Eqs . (|10p-( | l2|) scale as T /5. To maintain 
the validity of Eqs. (|10|)-(|l^y7however, one must impose 
the limit 5 < l\Eis, where A^;^ « 3x10^^ s"^ is the 



fine structure splitting between Pj 



1/2 



and ^P-iji- Physi- 



cally, this means that the light shifts due to the excited 
^^3/2 states remain negligible compared to those origi- 
nating from the states. As a concrete example, we 



consider the case b = 3x10 s~ and /q = 10 /g, where 
Is is the saturation intensity. This results in a pulse se- 
quence length of T « 35 /zs and a scattering time of an 
individual atom of T^c = l/Psc ~ 25T. We also men- 
tion that the ratio Tsc/T can be enhanced by pairing ev- 
ery beam in the original sequence with an appropriately 
chosen second beam. If both beams have the same po- 
larization and are detuned equally above and below the 
^"5*1/2^^^3/2 transition, the light shifts from the ^P3/2 
states are approximately cancelled out, while those from 
the ^Pi/2 states are reinforced. This makes it possible 
to abandon the requirement 5 <C AE'fg and choose much 
larger detunings. 

Finally, we note that the solution (^ for the dynamics 
of the system under the laser beams is only approximate 
and overlooks subtler but important points. Clearly such 
a transformation is unphysical in our case at the origin as 
it would introduce an infinite phase gradient energy 
and a more careful analysis of the dynamics near the 
center would have to include other terms in the Hamil- 
tonian. Also, in a real system, an absorption plate can- 
not perfectly generate the intensity profiles given in ( p^ ) 
and (16), which leads to various degrees of imperfection 
in vortex creation In addition, experimental efforts to 
create a soliton in a scalar condensate using an intensity 
profile like that of ( p7| ) show that the dipole potential also 
creates a density wave in addition to the soliton 
We refer the reader to the cited references for further 
discussion. Ultimately, however, it does not appear that 
these imperfections preclude the creation of vortices and 
solitons, which suggests that the creation of monopoles 
via phase imprinting should also be viable in spite of such 
imperfections. 



V. SUMMARY 

We have proposed a method to create a monopole in an 
antiferromagnetic spin-1 BEC by extending the phase im- 
printing concept used in scalar BEC's. We demonstrate 
that laser beams with different polarizations generate en- 
ergy shifts that can be used to engineer the phases of the 
magnetic substates in a broad manner. We also present 
an explicit pulse sequence that, starting from the ground 
state of the condensate, produces phases in the substates 
corresponding to those of a monopole. 
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Pulse 


Polarization 


Intensity profile 


Time (in units of n/ao) 


A 


a+ 


pi(r) 


1 


B 


a~ 


P2(r) 


1 


C 


a+ 


P3(r) 


1 


D 


TT 


P4(r) 


1/2 



TABLE I: The pulse sequence to create the monopole is shown, specifying the beam polarizations, intensity profiles, and pulse 
durations. Strictly speaking, the polarization and intensity profile of pulse C are inconsistent, but one can achieve the desired 
effect using a laser pulse inserted between two magnetic field pulses. 
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FIG. 1: The values for Cfj, scaled so that the smallest value is an integer. The solid lines give the coefficients for 7r-poIarize 
light, while the dashed lines correspond to o" "'"-polarized light. The values for (T~-polarized light are obtained from the a 
values by making the substitutions mp^ — mp in the diagram. 
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FIG. 2: The relevant frequencies of the atom-laser system. 



